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ABSTRACT 

A new method which enables us to estimate the inclination of Class I young stellar objects 
is proposed. Since Class I objects are not spherically symmetric, it is likely that the observed 
feature is sensitive to the inclination of the system. Thus, we construct a protostar model by 
carefully treating two-dimensional (2D) radiative transfer and radiative equilibrium. We show 
from the present 2D numerical simulations that the emergent luminosity £sedj which is the 
frequency integration of spectral energy distribution (SED), depends strongly on the inclination 
of the system i, whereas the peak flux is insensitive to i. Based on this result, we introduce a novel 
indicator /l, which is the ratio of -Lsed to the peak flux, as a good measure for the inclination. 
By using /l, we can determine the inclination regardless of the other physical parameters. The 
inclination would be determined by /l within the accuracy of ±5°, if the opening angle of bipolar 
outflows is specified by any other procedure. Since this spectrophotometric method is easier than 
a geometrical method or a full SED fitting method, this method could be a powerful tool to 
investigate the feature of protostars statistically with observational data which will be provided 
by future missions, such as SIRTF, ASTRO-F, and ALMA. 

Subject headings: ISM: clouds — radiative transfer — stars: formation — stars: pre-main-sequence 



1. Introduction 

Low mass pre-main-sequence stars are classified 
into four categories, Class 0, I, II, and III, mainly 
according to the shape of their spectral energy dis- 
tributions (SEDs) (Lada & Wilking 1984; Adams, 
Lada, & Shu 1987, Andre, Ward-Thompson, & 
Barsony 1993). The shape of SED is strongly de- 
pendent on the structure surrounding a central 
star, such as an envelope, an outflow, a disk, and 
so on. For example, since Class and Class I ob- 
jects show a peak in their SEDs in a sub- millimeter 
range, they are thought to be surrounded by a 
cool massive envelope, which heavily obscures a 
hot central part of the object. It is considered 
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generally that the order of the classification, 0, I, 
II, and III corresponds to an evolutional sequence 
of a low mass pre-main-sequence star. Class and 
I objects are thought to be protostars, which are 
enshrouded by a thick envelope, and to be younger 
than T Tauri stars, which are classified as Class II 
and III objects. 

A protostar is not a spherically symmetric sys- 
tem as implied by the presence of a bipolar out- 
flow. This asymmetry could come from a struc- 
ture such as a central star with a disk system, 
which seems to be the progenitor of a planetary 
system. Thus, a careful scrutiny of the struc- 
ture of a protostar system is important to under- 
stand the formation process of a star and/or plan- 
etary system. However, it is not easy to reveal 
the structure directly by imaging observations, be- 
cause optical imaging is significantly subject to the 
optically thick envelope and also radio imaging 
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is strongly limited by the resolution of the tech- 
nique. In contrast, a spectrophotometric method 
using SEDs seems to be effective to investigate a 
large number of protostar systems (Kenyon, Cal- 
vct, & Hartmann 1993; Men'shchikov, Hcnning, 
& Fischer 1999; Nakamoto & Kikuchi 1999). For 
instance, Kenyon et al. (1993) modeled a proto- 
star with incorporating spherically symmetric ra- 
diative transfer and simulated its SED. They fit- 
ted simulated SEDs with observed SEDs and in- 
ferred the structure of Class I objects. The accu- 
racy of the spectrophotometric method depends 
on the accuracy of the model of protostar, which 
is characterized by the density and temperature 
distributions in the envelope, the inclination of the 
system, and so forth. In particular, the inclination 
of the system is an important piece of information 
to probe the inner structure of the protostar. It 
can affect the SED significantly through 2D be- 
havior of radiative transfer. Hence, the inclina- 
tion is of essential importance in a realistic model 
of protostars. In other words, there may be a pos- 
sibility to estimate the inclination using the SED. 
So far, the inclination of a protostar has been es- 
timated by the geometrical feature of molecular 
outflow lobes. However, the geometrical method 
requires high spatial resolution of observations and 
long integration time. Thus, it is not effective for 
a statistical study. On the other hand, if we can 
estimate the inclination of the protostar solely by 
the SED, it would be quite effective to investigate 
the protostar structure. 

In this paper, we propose a new method to de- 
termine the inclination of Class I objects using 
their SEDs. We solve the radiative equilibrium 
for the 2D axisymmetric density distribution with 
incorporating 2D radiative transfer and simulate 
SEDs for the 2D protostar model. Then, it is 
shown that the ratio of the quantity £sed, which 
is the integration of the SED over frequencies, to 
the peak flux of the SED, (vL v ) ma x, can be a good 
measure of the inclination. The accuracy and the 
applicability of the present method is discussed as 
well. 

This paper is organized as follows. In section 
2 and 3, our model of a protostar and calculation 
procedure of the new method will be described. 
The new method will be introduced in section 4. 
The validity, reliability, and accuracy of the new 
method is also discussed in this section. The con- 



clusions of this paper will be given in section 5. 

2. Model 

The present two-dimensional axisymmetric 
model for a Class I object consists of three com- 
ponents; a central star, a circumstellar disk, and 
an envelope surrounding the former two compo- 
nents. The main opacity source is dust grains. It 
is assumed that the dust grains are well mixed 
with hydrogen molecule gas and the mass ratio 
to the hydrogen gas is 0.01. The opacity of dust 
grains as a function of the frequency is taken from 
Miyake & Nakagawa (1993). Physical parameters 
that characterize our model are summarized in 
Table 1. 

The central star is only the energy source and 
radiative heating is considered. Compressional 
heating, viscous heating, chemical heating, and 
any kind of heating sources other than the radia- 
tion from the central star are ignored, since they 
are negligibly small in almost all the regions of in- 
terest. The central star is assumed to have the 
surface temperature T* and emit the blackbody 
radiation with that temperature. The luminosity 
of the central star is given by L*. The mass of the 
central star is denoted by M*. 

Our model contains a circumstellar disk, al- 
though the presence of the disk has not been con- 
firmed observationally in Class I objects. Its in- 
ner and outer radii are assumed to be 0.1 AU and 
100AU, respectively. The surface density distribu- 
tion of the disk is given by E(r) = £i(r/lAU)~ p , 
where r is the radius in cylindrical coordinates 
from the central star, Si is surface density at 
1AU, and p is the index of the power, respectively. 
We adopt p = 1.5 which is a standard value for 
the minimum mass solar nebula model (Hayashi, 
Nakazawa, & Nakagawa 1985). The value of p 
does not affect the SEDs significantly. Vertical 
structure in the disk is determined by the hydro- 
static equilibrium between gravity by the central 
star and thermal pressure. 

We assume the outer radius of the envelope to 
be 1, 000 AU. Contributions to SED from the out- 
side of that region are neglected. This assumption 
seems appropriate, since the observed region of a 
Class I object is as small as a few thousand AU. 
In the envelope, the density distribution is given 
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by 

P {R) = pi (y^fj) . for 1AU < R < 1. OOOAU, 

(1) 

where i? is distance from the central star in spher- 
ical coordinates, p\ is mass density at 1AU, and q 
is the index of the power, respectively. Here, we 
adopt q = 1.5 because it is the value expected for 
the free-falling envelope (Shu 1977). Observation- 
ally, it is known that Class I objects usually have 
bipolar outflows, and it is suggested that the out- 
flow excavates the envelope and affects the SED 
considerably. In order to take into account the 
effects of the bipolar outflow, we situate a cavity 
along the symmetric axis. The cavity is assumed 
to be a conical shape with a half opening angle #bp, 
which is the angle between the symmetric axis and 
the boundary of the cavity. The density in the cav- 
ity is assumed to be 0.01 times smaller than the 
value which is given by eq.(l). Thus, the density 
in the envelope is given by 

P(W) = IS f 1 ^ 1 " (2) 

' K ' ' \ 0.01 x p(R) for 9 < 6» bp . v ' 

3. 2D Radiative Equilibrium Calculations 
3.1. Numerical Method 

In the first step, we obtain the temperature 
distribution in the model under the condition 
of the radiative equilibrium; X^ s J v dv = 
Jo°° X^Bvdv, where x^ bs is the absorption co- 
efficient, J v is the mean intensity, and B v is the 
Planck function, respectively. To calculate the ra- 
diative equilibrium, we solve 2D radiative trans- 
fer by a numerical scheme based on the so-called 
Variable Eddington Factor method (Stone, Mi- 
halas, & Norman 1992; Kikuchi, Nakamoto, & 
Ogochi 2002). The primary difference of our algo- 
rithm from that by Stone et al. is that we treat a 
frequency dependence of radiative transfer accu- 
rately. In our method, the radiation field is sepa- 
rated into two parts; one is the direct component 
from the central star and the other is the diffuse 
component. The direct component is calculated 
straightforwardly. The diffuse component is ob- 
tained by solving moment equations (Stone et al. 
1992), which are extended to a non-gray version 
in our method, coupled with the energy equation. 
To close the moment equations, the variable Ed- 



dington factor is introduced, which is determined 
by solving the radiative transfer equation taking 
the scattering effect into account. We expand the 
scattering phase function by the Legendre func- 
tion to the 2nd order. The set of equations are 
integrated with time. When the temporal change 
of the temperature distribution in the calculation 
region becomes small enough, we regard it as an 
equilibrium state. In the second step, the SED of 
the object is simulated by the direct ray-tracing 
with the radiative equilibrium temperature distri- 
bution. 

Our numerical code is based on cylindrical co- 
ordinates. We use 50 x 50 grids for r-z coordi- 
nates in space, and 50 x 100 grids for azimuthal 
and zenith angles to express directions in which 
the intensity propagates. Calculations with dou- 
ble number of grids in each dimension (16 times 
larger calculations) showed that the difference of 
results was within about 10% in /l, which will be 
defined in cq.(4). The frequency space is divided 
into 51 grids. 

3.2. SED calculation 

Simulated SEDs arc illustrated in Fig. 1. It 
is seen that SEDs are sensitive to i, particularly 
in the frequency range higher than the peak fre- 
quency fpeak- This change reflects the asymmetry 
of the density and temperature distribution. The 
optical depth along the line-of-sight from an ob- 
server to the central star, r c , is ~ 1 at ^ poa k, while 
r c > 1 for v > fpeak and t c < 1 for v < ^ p0 ak- 
Thus, for v < ^p 0a k, the dependence on the incli- 
nation is within a factor of 2 ~ 3 and the emergent 
flux is the superposition of the thermal radiation 
from dust grains, which is determined by the tem- 
perature and the total amount of dust, regardless 
of density distributions. In contrast, for v > t'peak, 
the anisotropy of flux is very large (2-4 orders of 
magnitude), and the emergent flux is dominated 
by the scattered radiation originating from the 
central star and inner disk, which is strongly af- 
fected by the density distribution of the envelope. 
Therefore, the flux in the pole-on (i = 0°) direc- 
tion is large, while the flux in the edge-on (i = 90°) 
direction is small. 
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3.3. Comparison between Semi and Full 
2D Calculations 

To see the importance of full 2D radiative trans- 
fer calculations, we performed approximate semi 
2D calculations, which were carried out by Kenyon 
ct al. (1993). They derived the temperature dis- 
tribution in the envelope from radiative equilib- 
rium, but they assumed the spherically symmetric 
density distribution when they calculated the ra- 
diative equilibrium, so they obtained spherically 
symmetric temperature distribution. To do that, 
first the spherically averaged density profile from 
the 2D axially symmetric density distribution is 
obtained, and then the ID spherically symmet- 
ric temperature profile is calculated. Using the 
spherical temperature distribution and the non- 
spherical 2D density distribution, the emergent 
SEDs from the object are obtained by ray-tracing 
. Following their procedure, we also calculated the 
SED from the semi 2D calculations and compared 
the results with our full 2D calculations. 

In Fig. 2, we compare two temperature distri- 
butions derived from two procedures. The temper- 
ature distribution obtained by the full 2D calcula- 
tions is not spherically symmetric at all. The most 
noticeable difference between two temperatures is 
found in the disk and the outflow regions. The 
temperature in the disk in the full 2D calculation 
is much lower than that of the spherical tempera- 
ture, because the disk has higher density than the 
spherically averaged density. We also find that the 
temperature in a shaded region, which is behind 
the disk with respect to the central star, is lower 
than the spherical temperature, though the den- 
sity in that region is not so different from spheri- 
cally averaged one. This is because that the total 
amount of radiation from the central star to the 
region is reduced by the efficient disk absorption. 
In contrast, the temperature in the outflow region 
in the full 2D calculation is higher than spheri- 
cal one, because the density in the region is lower 
than the spherically averaged one. Consequently, 
it is obvious that the temperature distribution ob- 
tained by the full 2D radiative equilibrium calcu- 
lation is quite different from spherically symmetric 
distribution. 

Emergent SEDs based on the temperature dis- 
tributions displayed in Fig. 2 are shown in Fig. 3. 
For each inclination, it is seen that (1) the flux de- 



rived from the full 2D calculations in a frequency 
range from 10 12 Hz to 10 13 Hz is lower than that of 
the semi 2D calculations, and (2) in a range from 
10 13 Hz to 10 15 Hz, the full 2D calculation flux is 
higher than the semi 2D calculation one. The first 
difference is due to the decrease of the envelope 
temperature compared to the spherical tempera- 
ture. Since the optical depth for the radiation at 
this frequency from the central star to the observer 
is around unity, the flux in this range is expected 
to be in proportion to the dust temperature in the 
envelope. On the other hand, the second difference 
is caused by the different mean intensity J„ espe- 
cially in the outflow region, because flux in this 
frequency range is mainly determined by the scat- 
tering, which is evaluated by y^^Jv, where xt° a 
is scattering coefficient, and is very anisotropic in 
the non-spherical density distribution cases. We 
can see from Figs. 2 and 3 that the full 2D radia- 
tive equilibrium calculations are indispensable to 
study the detailed structure of Class 0/1 objects 
using their SEDs and/or images. 

4. Estimation of Inclination Angle 
4.1. A New Inclination Indicator: /l 

The inclination dependence of SED provides a 
tool to infer the inclination itself. We introduce 
the emergent luminosity Lsed, which is defined 

by 

^SED = AttD 2 / {vF v )d\ogv, (3) 
Jo 

as a tracer for the change of SED, where D is 
the distance to the target object and F v is the 
observed flux at frequency v. Also, using the peak 
flux {vL v ) mayL = (v ■ AirD 2 F v ) max (flux at ^ poak ), 
we define a ratio /l by 

, _ £SED 

(Note that Lsed, {vL v ) max , and /l arc all ob- 
servable quantities. We can obtain them from ob- 
served SEDs.) We show the change of /l against 
i in Fig. 4. Fig. 4(a) indicates that i can be 
determined if /l is evaluated from observed SED 
for the angle of i > 9b P - But, when i is in a range 
i < #b P , /l is almost independent of i. This means 
that SED does not change if we observe a proto- 
star through the bipolar cavity. For i ~ 90°, /l is 
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also independent of i. It corresponds to the case 
that we observe a protostar through the disk. 

The level of /l depends on 9b P - It implies that 
$b P is one of the important parameters to deter- 
mine the anisotropy of a protostar system. The 
ratio /l increases for all i as #b P increases, and 
it is worth noting that there exist maximum and 
minimum of /l for each #bp- 

We calculated /l for full 2D and semi 2D results 
shown in Figs. 2 and 3 to clarify the importance of 
the full 2D calculation. These results are shown in 
Fig. 5. In the semi 2D calculation, the value of /l 
is lower than that of the full 2D calculation, espe- 
cially for the small inclination. It is because that 
the peak fluxes (^L„) max are different by a factor 
of ~ 2 while the difference of observed luminosi- 
ties Lsed is so small. For the large inclination, 
the differences of (^L„) max and Lsed are canceled 
each other so that the values of /l are almost the 
same. 

4.2. Test of the Method 

We test our new method with a real object 
TMC1A (IRAS 04365+2535). TMC1A is a Class 
I object (e.g. Kenyon et al. 1993; Chandler et al. 
1998) and has an outflow with the half opening 
angle 6> bp ~ 20° (Chandler et al. 1996; Gomez 
et al. 1997). According to Motte et al. (2001), 
TMC1A has an extended envelope and is recog- 
nized as a 'true protostar' not a 'transition ob- 
ject' nor an 'edge-on view of T Tauri star'. Thus, 
we think that TMC1A is an archetype of Class 
I objects and an appropriate example to test our 
new method to assess the inclination. Observa- 
tional data of TMC1A are taken from Myers et al. 
(1987), Kenyon et al. (1993), and Chandler et al. 
(1998). 

Here, we adopt the new method to derive the 
inclination, and compare the estimate with the re- 
sults by the full SED fitting. First, we attempt to 
estimate the inclination of the object using our 
new method. We evaluate first the observed lu- 
minosity Lsed by integrating the SED and obtain 
Lsed = 2.8Z/0, which is consistent with the results 
obtained by Myers et al. (1987) (L = 2.4L Q ) and 
Kenyon & Hartmann (1995) (L = 2.2L Q ). We 
next calculate /l from Lsed and (^L„) max , and 
obtain /l = 2.52. From Fig. 4(a), this value in- 
dicates that inclination angle is about 20°. This 



inclination is estimated with our new method only 
using the Lsed and the peak flux, without fitting 
all the physical parameters of the object. 

Next, we try to obtain all the physical parame- 
ters of the object by fitting the SED and estimate 
the inclination angle. In practice, we fit the SED 
going through the following steps: 

1. Set the outflow half opening angle 8b P = 20°. 

2. Infer the luminosity of the central star L* 
from the peak flux of the SED. 

3. Infer the circumstcllar mass M onv + ilLjisk 
(M env oc pi, Mdisk oc Ei) from the flux in a 
frequency range lower than the frequency of 
the peak flux f poa k. 

4. Infer the ratio of the envelope mass to the 
disk mass and inclination from the shape of 
the SED in a frequency range higher than 

^peak • 

The best fitted SED is displayed in Fig. 6 and 
adopted parameters are listed in Table 2. It is seen 
that our model can reproduce the observational 
data very well. In this procedure, the inclination 
of the object is estimated to be 22°. 

It is clear that the estimated values obtained by 
above two different methods agree quite well with 
each other. Since the estimation based on the full 
SED fitting is considered to be more reliable, it 
seems that the new spectrophotometric method 
provides an effective tool to assess the inclination. 

However, it should be noted that the inclina- 
tion angle obtained above does not agree with 
the values obtained by Kenyon et al. (1993), 
which is i ~ 60°, and Chandler et al. (1996), 
i = 40° — 70°. Such discrepancy could be at- 
tributed to the method used to infer the inclina- 
tion. Kenyon et al. (1993) used the spherically 
symmetric temperature distribution when they es- 
timated the inclination with the SED. But the 
spherical temperature distribution can lead a sig- 
nificant difference in the resultant SED as shown 
in section 3.3. Chandler et al. (1996) estimated 
the inclination angle based on the shape of outflow 
lobe on the sky. They first restricted the inclina- 
tion angle to be less than 70°. Then, they placed 
the lower limit of the inclination to be 40°, though 
the reason why they adopted the value was not 
given clearly. 
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4.3. Effects of Other Parameters 

Here, we estimate uncertainties of the present 
method by changing parameters of the model. As 
shown in Fig. 4(a), an uncertainty of /l is large 
if the half opening angle #b P is thoroughly un- 
known. But, it is possible to pose a constraint 
for 8 hp in the Class I phase by the observation 
of outflow lobes (Cabrit & Bertout 1986). For 
example, it is estimated that 6*b p = 13° - 19° 
for Class I object TMC1 (IRAS04381+2540) and 
15° - 21° for TMC1A (IRAS04365+2535), respec- 
tively (Chandler ct al. 1996). There arc also some 
other results for estimation of 6*b p : #b P = 22.5° for 
L1448 (Bachiller ct al. 1995) and B335 (Hirano 
ct al. 1988; Chandler & Sargent 1993), and 35° 
for L1551 IRS5 (Andre et al. 1990). From these 
estimations, it seems possible to constrain #b p to 
be 10° - 40°. This constraint would be justified by 
the fact that if #b P becomes larger than this range, 
resultant SED exhibits double peaked feature in 
far-IR and near-IR to optical, but these SEDs are 
not classified Class I category. From Fig. 4(a), a 
typical error for determination of i is ±10° if we 
restrict ourselves to #b P = 10° - 40°. 

Next, we examine to what extent the luminos- 
ity of the central star L» affects the ratio /l. Fig. 
7 shows calculated /l as functions of the incli- 
nation i with three different luminosities, O.ILq, 
ILq, and 1QL & . It is seen that when the luminos- 
ity of the central star is low (L* = 0.1i Q ), /l for 
the small inclinations decreases compared to the 
standard case (£* = ILq), but /l for the large 
inclinations does not change. In contrast, when 
the luminosity is high (L» = 10i Q ), /l for the 
large inclinations becomes higher than the stan- 
dard case, while /l for the low inclinations remains 
the standard value. It is true that /l is affected by 
the central star luminosity to some extent, but the 
change of the value in a range sensitive to the in- 
clination is not large. Thus, when we estimate the 
inclination of an object using this new method and 
/l, the effect of the central star luminosity seems 
to be negligible. 

Also, pi could affect the results. Fig. 8 shows 
the condition of /l = constant for each p\ in the 
i — @bp plane. The flux at v > i; poa k decreases 
with increase of pi, because r c is proportional to 
pi and the flux in this frequency range is pro- 
portional to exp(— t c ). In contrast, the flux at 



v < f P eak increases with increase of p\ because the 
total amount of dust grains is increased. Thus, /l 
tends to decrease with increase of p\ , so that the 
line of /l = constant is shifted above with increas- 
ing pi in the i~0b p plane. For 9 bp = 10° -40°, the 
inclination i changes by ~ 10° when pi changes by 
an order. 

By contrast, Si does not influence /l- Fig. 9 
displays calculated /l with different Si. It is eas- 
ily seen that Si does not change /l significantly. 
This is because the circumstellar disk contribute 
little to the column density along the line of sight, 
on which /l strongly depends. 

Therefore, we conclude that the primary un- 
certainty of our method comes from the informa- 
tion of $b P - If &bp can be estimated by any other 
means, then the error of inclination obtained by 
our method becomes roughly ±5°. If #b P is un- 
known, then the error may be roughly ±15°. 

5. Conclusions 

Using 2D radiation transfer calculations, we 
have obtained the radiation fields for 2D axisym- 
metric protostar model, consistent with the cen- 
tral star, the circumstellar disk, and the envelope, 
and derived the SEDs of Class I objects. We have 
found that the ratio /l between the emergent lu- 
minosity, isED, and the peak flux in the SED, 
(^i J/ ) max , is a good indicator of the inclination 
angle of the object, because Lsed is sensitive to 
the inclination, while the peak flux is insensitive. 
Through the test with real data and the analy- 
sis on the other physical parameters, it has been 
shown that /l is a robust tool to assess the incli- 
nation of a Class I object. Hence, /l can provide 
a new spcctrophotomctric method to estimate the 
inclination angle of a Class I object. It is benefi- 
cial that, in this new method, only the luminosity 
Lsed and the peak flux {vL v ) mayi are needed to 
estimate the inclination. The typical error of the 
method is roughly ±5° if the half opening angle 
#b P is known and roughly ±15° if it is unknown. 

The present method is applicable for a great 
deal of data provided by future missions of in- 
frared, sub-millimeter, and millimeter observato- 
ries, e.g., SIRTF, ASTRO-F, and ALMA. These 
missions are expected to reveal numerous pro- 
tostar candidates. Then, it may be possible to 
carry out precise statistical study of protostars. 



6 



The present spectrophotometric method could be 
a powerful tool in the statistical study. 
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in Center for Computational Physics, University 
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Symbol 


Meanings 


Parameter Range 


L* 


luminosity of central star 


0.1 - 1OL 


T* 


temperature of central star 


4000K 


M, 


mass of central star 


O.5M 


Si 


surface density of disk at 1AU 


200 - 15000 g cnT 2 a 


P 


power law index of disk 


1.5 


Pi 


density of envelope at 1AU 


10 -13.5_ 10 -12.5 gcm -3b 


q 


power law index of envelope 


1.5 




opening angle of bipolar cavity 


0°-55° 


i 


inclination angle 


0°-90° 



Table 1: Parameters for the present protostar 
model. 



"This corresponds to the disk total mass ~ 0.03Mq assum- 
ing dust-to-gas mass ratio of 0.01. 

b pi = 10~ 13 g cm~ 3 corresponds to the envelope total mass 
~ O.O3M0 assuming dust-to-gas mass ratio of 0.01 and 

e bp = o° 
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Object 


L* [L Q ] 


Si [g cm 2 ] 


Pi [g cm 3 ] 




i 


TO 


KCH 


Outflow 


TMC1A 


1.0 


10 4 


10 -ia.b 


20° 


22° 


20° 


60° 


40° - 70° 



Table 2: Best fit parameters for TMC1A (IRAS 
04365+2535). From second to fifth columns show 
best fit parameters in our model for TMC1A. 
The sixth column shows the inclination angle es- 
timated from /l. Seventh and eighth columns 
show the inclinations estimated from other work 
by Kenyon ct al. (1993; marked 'KCH'), who car- 
ried out full SED fitting using semi 2D protostar 
model, and Chandler et al. (1996; marked 'Out- 
flow'), who infered the inclination from the shape 
of outflow lobes on the sky. 
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Fig. 1. — Emergent spectral energy distribution 
(SED) with the present protostar model. Pro- 
tostar model parameters are L„ = IL Q , p\ = 
1(T 13 g cm- 3 , Si = 2000 g cm" 2 , and 9 hp = 25°. 
Four curves represent SEDs with the inclination 
i = 0° (solid curve), 30° (dashed curve), 60° 
(dotted curve), and 90° (dot-dashed curve), respec- 
tively. It is seen that the SED changes with the 
inclination i, even though the protostar structure 
does not change. The emergent luminosity £sed 
also changes with the inclination. 



Fig. 2. — Temperature distributions obtained by 
the full 2D radiative equilibrium calculation (solid 
curve) and the spherically averaged ID calcula- 
tion (dotted curve) in regions (a) 30 AU x 30 
AU and (b) 1000 AU x 1000 AU. Adopted pa- 
rameters are £* = 1.0L Q , pi = lO -13 g cm- 3 , 
Si = 2000 g cm -2 , and 9 hp = 25°. In the spheri- 
cally averaged ID calculation, density distribution 
is spherically averaged first, and then, the radia- 
tive equilibrium is solved with the ID spherically 
symmetric radiative transfer calculation. 
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Fig. 3. — Spectral energy distributions based on 
temperature distributions obtained by two differ- 
ent methods shown in Fig. 2 for different incli- 
nation angles. SEDs from the non-spherical tem- 
perature distribution (obtained by the full 2D ra- 
diative equilibrium calculation) are shown by solid 
curves, and those from the spherical temperature 
are shown by dotted curves. 
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Fig. 4. — Three observables obtained directly from 
SED are plotted against the inclination. In each 
panel, two different curves correspond to SEDs 
with different opening angles in the outflow region, 
6» bp = 20° (solid curve), and 40° (dotted curve), 
respectively: (a) the observable inclination indica- 
tor /l, which is the ratio of Lsed to {vL„) max , (b) 
the emergent luminosity £sed, an( ^ ( c ) ^ ne P ea k 
flux (i/L v ) max . It is seen that /l changes abruptly 
around i = 6^ p . If /l is obtained from the SED, 
the inclination i can be estimated for i > #bp by 
using this figure. 
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Fig. 5. — The inclination indicator /l based on the 
temperature distributions obtained by the semi 2D 
(dashed line) and the full 2D calculation (solid 
line). The values of /l are different from each 
other especially for the small inclination. 




Fig. 6.— Best fitted SED with the full SED fit- 
ting procedure for TMC1A (IRAS 04365+2535). 
Adopted parameters are as follows: = 1^0, 
pi = 10~ 12 ' 5 g cm~ 3 , Ei = 10 4 g cm~ 2 , and 
#b P = 20°, and i — 22°. Observational data of 
the object, plotted by filled triangles, are taken 
from Myers et al. (1987), Kenyon et al. (1993), 
and Chandler et al. (1998). 



0. 1 Lsun 
1 Lsun 
10 Lsun 




10 20 30 40 50 

j [deg] 



70 80 90 



Fig. 7. — The inclination indicator, /l, vs the 
inclination, i, with three different luminosities. 
When the luminosity is low, /l for the small incli- 
nation is reduced, while if the luminosity is high, 
/l for large inclination is raised. 
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Fig. 8. — The condition of /l = 3 is plotted on i- 
f?bp plane. Three curves correspond to three den- 
sities in the envelope at 1AU, pi = 10~ 12 ' 5 g cm~ 3 
(solid curve), pi = 10~ 13 -°g cm~ 3 (dashed curve), 
and pi — 10~ 13 5 g cm -3 (dotted curve), respec- 
tively. An error due to the change of pi is esti- 
mated to ±5° in the range of 6» bp = 0° - 40°. 
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Fig. 9. — The inclination indicator, /l, vs the in- 
clination, i, with three different disk masses. The 
disk mass does not affect the ratio /l significantly. 
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